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Summation of the Virial Expansion for Adsorption Isotherm in the Cayley Tree 
Approximation 


WSTĘP 


Celem tej pracy będą badania teoretyczne adsorpcji gazów na ciałach 
stałych, oparte na nowoczesnym formalizmie opisu tego zjwiska, polega- 
jącym na rozwijaniu odpowiednich wielkości termodynamicznych w sze- 
regi według potęg małego parametru. Najbardziej znane są standartowe 
rozwinięcia w potęgi aktywności (aktywnej gęstości), stanowiące podsta- 
wową formę rozwinięcia. Za formę pochodną można uważać rozwinięcia 
wirialne rozważane zgodnie z tradycją za rozwinięcia w potęgi gęstości. 
Szeregi wirialne bowiem uzyskuje się drogą odpowiednich przekształceń 
szeregów standartowych. Użyteczność takich rozwinięć ogranicza się oczy- 
wiście do układów, w których zarówno aktywność, jak i gęstość stanowią 
rzeczywiście małe parametry. Są to więc przede wszystkim układy gazo- 
we o niskiej gęstości i umiarkowanej temperaturze. Wówczas bowiem 
można zadowolić się kilkoma pierwszymi wyrazami szeregu, uważając że 
wyrazy z wyższymi potęgami wnoszą niewielki wkład do własności ma- 
kroskopowych układu. 

Oprócz wspomnianych wyżej rozwinięć standartowych istnieją inne 
rozwinięcia, z których najbardziej znane jest rozwinięcie Urselła-Mayera. 
Rozwinięcie to w najbardziej pierwotnej postaci jest rozwinięciem ciśnie- 
nia i średniej liczby cząstek w potęgi aktywności. Można je zatem rozpa- 
trywać również jako metodę analizy współczynników rozwinięć standar- 
towych. Ta metoda analizy, którą będziemy dalej nazywać metodą Ur- 
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sella-Mayera, stosowana była najpierw do układów monocząstkowych, 
w których energia potencjalna pochodziła wyłącznie z oddziaływań mię- 
dzycząstkowych. 

Dopiero w latach pięćdziesiątych pojawiły się próby jej adaptacji do 
układów monocząstkowych, znajdujących się w zewnętrznym polu po- 
tencjalnym. Rzecz ciekawa, próby te pojawiły się właśnie w związku z 
opisem teoretycznym faz adsorpcyjnych. 

Pionierem pod tym względem był Ono. [2, 3], autor licznych prac, 
które zapoczątkowały prawdziwie nowoczesny opis faz powierzchniowych 
adsorpcyjnych. W pracach tych traktuje się fazę adsorpcyjną jako pewien 
gaz rzeczywisty w polu zewnętrznym sił adsorbentu. W dalszych pracach 
Ono [4] uogólnia swoje rozważania na przypadek układów heterocząst- 
kowych. i 

W nastepnych latach ukazuja sie prace rozwijajace ten punkt widze- 
nia na zjawisko adsorpcji z gazów [5—11]. Byly to w przewazajacej mie- 
rze prace oparte na formalizmie termodynamiki statystycznej klasycznej, 
czego usprawiedliwieniem jest fakt, że większość układów adsorpcyjnych 
o znaczeniu praktycznym to układy z cząstkami adsorptywu stosują- 
cymi się dobrze do praw mechaniki klasycznej. Pojawiły się też prace o 
adsorpcji gazów „kwantowych” [12]. Zasadniczą cechą wspólną większoś- 
ci tych publikacji była procedura obcinania rozwinięć na określonej potę- 
dze małego parametru. Była to zresztą naturalna konsekwencja stosowa- 
nia podobnej procedury w opisie zwykłych układów gazowych, bez obec- 
ności pola zewnętrznego. 

Tymczasem w opisie układów monocząstkowych bez pola zewnę- 
trznego pojawiła się nowa idea, której autorem jest Fuliński. Idea ta 
zasadza się na sumowaniu, zarówno rozwinięć standartowych [13], jak 
i innego typu [14] aż do nieskończoności. Oczywiście zupełnie dokładne 
sumowanie tego typu jest niemożliwe. Z każdego wyrazu rozwinięcia wy- 
biera się tylko pewien wkład. W metodzie Fulińskiego można wyróżnić 
trzy zasadnicze etapy: 

1) analizę współczynników rozwinięć, bądź to w oparciu o dotychcza- 
sowe metody, bądź w oparciu o zupełnie nowe własne metody: 

2) rozbicie wartości współczynników na pewne klasy wartości i ocenę 
wkładu tych klas; 

3) podanie metody sumowania rozwinięcia, przy założeniu wkładu od 
pewnych tylko klas. 

Postępowanie powyższe okazało się bardzo owocne. Uzyskane wyniki 
zawierały przepowiednie istnienia przejść fazowych pierwszego i drugie- 
go rodzaju. Autorów tej pracy od dłuższego czasu intrygował problem, czy 
idei częściowych sumowań aż do nieskończoności nie można by wykorzy- 
stać na gruncie termodynamiki statystycznej adsorpcji z gazów. Oczywiś- 
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cie zmieni sie tu stopien komplikacji oraz cele zainteresowan. Fazy ad- 
sorpcyjne trzeba rozważać jako znajdujące się w niejednorodnym polu 
zewnętrznym adsorbentu, zaś wielkością termodynamiczną najbardziej 
interesującą będzie tu nie ciśnienie, lecz średnia liczba czestek w ukła- 
dzie, ściślej biorąc — w fazie adsorpcyjnej. Ta średnia liczba cząstek ma 
oczywisty związek z izotermą adsorpcji — wielkością najczęściej badaną 
w praktyce. 

W porównaniu z fazą objętościową bez obecności pola zewnętrznego, 
rośnie tu bardzo stopień komplikacji, choć pojawiają się też fakty pozy- 
tywne z tego punktu widzenia. Do takich należy przede wszystkim fakt 
krótkiego zasięgu sił adsorpcyjnych i ich stosunkowo duża wartość. Te po- 
zytywne fakty stworzyły możliwość pewnych uproszczeń oraz możliwość 
wykorzystania pewnych wyników, osiągniętych wcześniej przez Fuliń- 
skiego. W niniejszej pracy zaprezentowano próbę częściowego sumowania 
rozwinięcia standartowego dla izotermy adsorpcji. Próba ta bazuje na me- 
todzie Ursella-Mayera, a celem jej jest przebadanie sumy rozwinięcia 
przy założeniu wkładu klasy drzewek Cayleya. 


CZĘŚĆ TEORETYCZNA 


Rozważania części teoretycznej rozpoczniemy od zdefiniowania kilku 
pojęć, niezbędnych w dalszym wykładzie zagadnienia. Są to dobrze zna- 
ne pojęcia z teorii grafów liniowych Ursella-Mayera [15]. Należy zwró- 
cić uwagę, że przyjęta nomenklatura jest robocza z uwagi na dotychcza- 
sowy brak oficjalnych ustaleń w tej dziedzinie. i 

Grafem liniowym nazywamy skończony zbiór P punktów oraz skon- 
czony zbiór L linii, z których każda łączy dwa różne punkty zbioru P. 
Każda para punktów nie jest połączona więcej niż jedną linią. Jeśli za- 
równo P, jak i L są zbiorami pustymi, wówczas G nazywamy grafem pu- 
stym i oznaczamy go zazwyczaj przez Q. 

Jeśli G i H są takimi grafami, że H można otrzymać z G drogą wyję- 
cia ze zbioru G pewnego podzbioru jego punktów, wraz ze wszystkimi 
liniami wychodzącymi z jednego z punktów tego podzbioru, wówczas H 
nazywamy podgrafem grafu G. 

Jeśli dalej G i H są takimi grafami, że H można otrzymać drogą wy- 
jęcia jego linii, wówczas graf H nazywany jest grafem prostszym od G. 
Łańcuchem nazywamy sekwencję punktów (Po, Py, .., Pn), gdzie n = 1, 
taką, że dwa następujące po sobie punkty sekwencji spięte są linią. Dwa 
punkty P i Q spięte są za pomocą łańcucha, jeśli istnieje łańcuch (Po, Pi, 
PIES Q. 

Graf nazywamy spójnym, jeśli każda para jego punktów jest spięta za 
pomocą łańcucha. Graf zawierający jeden tylko punkt definiuje się jako 
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spójny, ale graf pusty definiuje sie jako niespójny. Spójnością x=x (G) 
grafu nazywamy minimalną liczbę punktów, których wyjęcie wraz z linia- 
mi z nich wychodzącymi daje w efekcie graf niespójny. Zgodnie z powyż- 
szą definicją grafy niespójne mają x=0, a grafy spójne mają x» 1. 

Graf nazywany jest podwójnie spójnym lub nieredukowalnym, jeśli 
» =2. Grafem minimalnie spójnym nazywamy graf, w którym wyjęcie 
jednej choćby linii daje w efekcie graf niespójny. Podobnie grafem mi- 
nimalnie podwójnie spójnym nazywamy graf, w którym wyjęcie której- 
kolwiek z jego lini prowadzi do nowego grafu, który nie jest już podwój- 
nie spójny. 

Drzewem Cayleya nazywamy taki graf, w którym każda para jego 
punktów spięta jest najwyżej jednym łańcuchem. Gwiazdą nazywamy 
pcdwójnie spójny graf, zawierający co najmniej trzy punkty i liczbę linii 
równą liczbie punktów. Drzewem Hushimi nazywamy taki graf spójny, 
w którym maksymalnie nieredukowalnym podgrafem jest gwiazda lub 
„pojedyncza linia”, tzn. graf składający się z dwu punktów oraz jednej 
linii. 

Jak już we wstępie stwierdziliśmy, metoda Ursella-Mayera jest sposo- 
bem analizy współczynników rozwinięć. Metoda ta pozwala podzielić każ- 
dy współczynnik na pewną sumę prostszych wyrażeń, będących na ogół 
calkami wielokrotnymi. Z każdą z całek związany jest pewien graf Maye- 
ra, charakteryzujący zupełnie postać całki, do której się odnosi. Całki te 
z powodu pewnych analogii nazywane są całkami konfiguracyjnymi, choć 
nie jest to zupełnie słuszne. Można to sformułowanie jednak przyjąć z za- 
strzeżeniem, że będzie używane dalej tylko w odniesieniu do metody Ur- 
sella-Mayera. ~ 

Wartością grafu nazywamy wartość całki konfiguracyjnej, do której 
się graf odnosi. W ten sposób można przejść od rozwirięcia według potęg 
do rozwinięcia według grafów. Przez rozwinięcie graficzne będziemy da- 
lej uważać szerez, którego wyrazy są wartościami pojedynczych grafów. 
Jeśli szereg taki jest zbieżny, wówczas mówi się, że jest on równy sumie 
wkładów od wszystkich grafów szeregu. 

Rozwinięcia graficzne można klasyfikować według ich stopnia d. Stop- 
niem graficznego rozwinięcia nazywamy minimalną spójność grafów wy- 
stępujących w rozwinięciu. W badaniach prowadzonych dotychczas roz- 
ważano graficzne rczwinięcia stopnia zerowego, pierwszego, drugiego 
i trzeciego. 

Wprowadźmy jeszcze następujące oznaczenia. Zbiór d-spójnych gra- 
fów, zawierających s punktów oznaczamy przez I gl», zaś zbiór wszystkich 
Ig, niezależnie od s oznaczamy przez l'a. Zwyczajowo przez Acr ozna- 
cza się zbiór drzew Cayleya, a przez Acr™ podzbiór składający się z drzew 
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Cayleya zawierających s punktów. Oznacza się dalej przez £a podzbiór 
wszystkich prostszych grafów ze zbioru Ig. Innymi słowy, Za składa sie 
ze wszystkich grafów G ze zbioru Ig takich, ze nie istnieje żaden inny 
graf H w zbiorze Ia, który jest prostszy niż G. Analogicznie przez Za 
oznaczamy zbiór wszystkich prostszych grafów w Iy. Z powyższych 
oznaczeń wynika, ze X,(0 — Ac4(), oraz Y= Acr. 

Pewien zbiór grafów, których suma wielkości jest równa pewnej wiel- 
kości A oznacza się jako I'a(4). Odpowiednio podzbiór wszystkich prost- 
szych grafów ze zbioru l'a (A) oznaczamy przez Xa (A). 

Jak powszechnie wiadomo, obliczanie i sumowanie wartości grafów ja- 
kiegokolwiek zbioru Ig (A) jest obecnie niemożliwe. Istnieje jednak spo- 
ro metod przybliżonego obliczania określonej wielkości fizycznej A. Moż- 
na je rozpatrywać jako procedury polegające na zsumowaniu pewnego 
podzbioru Dl'49(A) ze zbioru Ta (A). Procedura Fulińskiego polega na su- 
mowaniu wartości grafów pewnego zbioru Za (A), przy czym liczba ele- 
mentów zbioru jest nieskończona. Jest to istotna różnica w porównaniu 
z najczęściej dotychczas stosowanymi metodami aproksymacyjnymi, któ- 
re obcinały rozwinięcia na pewnych niskich potęgach małego parametru. 
Jest to równoważne z sumowaniem zbiorów o skończonej liczbie elemen- 
tów. Okazało się przy tym, że dopiero sumowanie zbiorów o nieskończo- 
nej liczbie grafów prowadzi do opisu tak interesujących zjawisk krytycz- 
nych. Wiąże się to, zdaniem Fulińskiego, z tym, że tego rodzaju sumowa- 
nie opisuje długozasięgowe oddziaływania, odpowiedzialne za istnienie 
zjawisk krytycznych. Pewne powodzenie teorii Percusa-Yevicka i „Hy- 
pernetted Chain” teorii wiąże się prawdopodobnie z tym, że są to proce- 
dury prowadzące w istocie do sumowania klas grafów o nieskończonej 
liczbie elementów. Niestety, ich związek z rozwinięciem graficznym Ur- 
sella-Mayera, a co z tym się wiąże, i stopień przybliżenia nie zostały do- 
tąd wyjaśnione. 

W niniejszej pracy dokonamy próby sumowania wartości zbioru 
Acr (N), gdzie N jest średnią liczbą zaadsorbowanych cząstek. Rozważmy 
w tym celu pewien układ adsorpcyjny o ogólnej objętości Vr. Układ ten 
możemy podzielić na trzy strefy objętościowe; Va, V, i V; zdefiniowane 
następująco: 

1) V, jest objętością zajmowaną przez masę adsorbentu; 

2) Vs jest częścią objętości zajmowanej przez adsorptyw, taką, że ist- 
nieją w niej znaczne siły adsorpcji; 

3) V; jest częścią objętości fazy adsorptywu, w której praktycznie nie 
istnieją siły adsorpcji. 

Oczywiście Vq= V44- Vs+ V;. 

Definiujemy obecnie potencjał adsorpcji s następująco: 
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+05 dla PFE K 
Er) = Edd Je iG (1) 
0 da TEX 


Powyżej r jest wektorem wcdzącym środka masy molekuły. Dla wiek- 
szości układów adsorpcyjnych spotykanych w praktyce można założyć z 
dobrym przybliżeniem plararność powierzchni międzyfazowej gaz-ciało 
_ stałe. Wówczas potencjał e ( r) można przedstawić w innej jeszcze postaci. 
Niech płaszczyzna xy będzie płaszczyzną adsorbentu. Wówczas 


*oo dia 0£z<a 
E(F)=€(z)=i6, dla agz<b (2) 
0 dia b«&z««*o 


Przyjmować dalej będziemy, ze potencjał oddziaływań międzycząstko- 
wych U (m, Ta) można przedstawić w postaci: 


+0o dla O«r <6 
UIR,RI=UIF)=4U dla 6<r<R6 (3) 
0 dla Ré&r«*o 


gdzie r=|r; — ral 

Oczywiście istnieją bardziej realistyczne postacie potencjału, zarówno dla 
e (z) [16], jak i dla U (r) [17]. Wiadomo jednak również, że potencjały stu- 
dniowe ze wzorów (2) i (3) są dobrym przybliżeniem przy odpowiednio do- 
branych parametrach. 

Rozważmy obecnie nowy hipotetyczny układ, analogiczny do poprzed- 
niego, z tym jednak, że brak w nim sił adsorpcji, co zapiszemy iS 
pująco: 

ago ae dla O«z«a (b 

0 dla a£z« *oo 
gdzie e*(z) jest potencjałem dla powyższego hipotetycznego układu. Obli- 
czmy obecnie ilość Na zaadsorbowanych cząstek. Stosując metodę Ursella- 
-Mayera, uogólnioną na przypadek pola zewnętrznego, otrzymujemy [18]: 
N= ) Ub-BJA' (5) 
(x1 

We wzorze powyższym b, i bi* są calkami gronowymi, identycznymi do 
tych w zwykłym rozwinięciu Ursella-Mayera, z wyjątkiem faktorów 


fools] s een] 


i^i 
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pod całkami. Całkowanie odbywa się w nich po całym obszarze fazy ad- 
sorptywu (Vs+ V;). Dalej 4 jest aktywnością (aktywną gęstością). War- 
tości funkcji podcalkowych w całkach (bi bi*) są różne od zera tylko 
wtedy, gdy choć jeden z wektorów wodzacych Tri względem których od- 
bywa się całkowanie, znajduje się w obszarze V,. Dla większości ukła- 
dów adsorpcyjnych spotykanych w praktyce faktor exp —sy/kT jest wię- 
kszy niż 105 w zwykłych temperaturach. Wynika stąd, że wartość funkcji 
podcałkowej w calce (b; b,*) maleje bardzo szybko, w miarę jak rośnie 
liczba r; € Vy. Z tego też powodu można z bardzo dobrym przybliżeniem 
napisać, że: 


Ne= ) lib B) (6) 
Uv , 
gdzie (bi,— bis*) jest całką, w której całkowanie odbywa sie po obszarze 
V,. Mamy zatem: 


N= ) ut [exe| ape (7) 


W powyższym wzorze bu? są zwykłymi całkami gronowymi Ursella-Ma- 
yera, w których całkowanie ograniczono do obszaru V;. Równanie (6) 
można napisać w innej jeszcze postaci: 


N=) LBA, (8) 
(94 

gdzie a= [exp (—e,/kT)—1] A jest zrenormalizowaną aktywnością. Równa- 
nie (8) jest formalnie analogiczne do odpowiedniego równania opisujacego 
zwykłą fazę objętościową. Dlatego możemy obecnie skorzystać z gotowe- 
go wyniku otrzymanego przez Fulińskiego [13]. Jeśli w całkach 
bis® uwzględnimy wkłady pochodzące tylko od drzew Cayleya, wówczas 
na średnią liczbę cząstek Na, zaadsorbowanych w obszarze Vs, otrzymamy 
następujący wynik: 


N= [atA [0,12.8,- A.A]! (9) 
gdzie funkcje p-m (x) maja między innymi następujące własności: 
2 
Q(x) =,lx)- 4 [9,09] -X (10) 


Pix) =x exp[g tx] s 


Wielkość f, zdefiniowana jest następująco: 


ATY Uer] nj non (12) 
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g = 1075. 10751077 
| -6= 40% 


WMC 
UNE 
+ 
o 


Ryc. 1 


Wynika stąd, ze wielkość f, jest proporcjonalna do drugiego współczyn: 
nika wiriału gaz-ciało stałe. Stosując własności funkcji q—4 (x) otrzy- 
mujemy: 
= Na expl_g, Ne 
A K exp) Bs Ne) (13) 
Powracając do zwykłej aktywności A i wprowadzając pojęcie gęstości w 
fazie powierzchniowej g;=N,/V; otrzymujemy: 
= -&) ? 
A [eo (e p. expe B; fr) (14) 
Aktywność X. fazy zaadsorbowanej musi być oczywiście równa aktywności 
wolnej fazy adsorptywu. Wprowadzając pojęcie gęstości ọọ wolnej fazy 
adsorptywu mamy równocześnie: 


A=f exp A) (15) 


Sumowanie rozwiniecia wirialnego dla izotermy adsorpcii... 75 


Objętość wolnej fazy adsorptywu V; jest zazwyczaj dużo większa niż V; 
i na tyle duża, że można ją uważać za zwykłą fazę objętościową. Wówczas: 


UIT) = 
A- | few A ]-1|a7 (16) 
v 
Porównujac równania (14) i (15) otrzymujemy: 


R [exo (£ )]exptp p= Po &xP(-Bo fà) (17) 


Równanie (17) można traktować jak nowe równanie izotermy adsorpcji. 
Wiąże ono bowiem gęstość w fazie powierzchniowej z gęstością wolnej 
fazy adsorptywu. Interesujący problem zbadania własności takiego rów- 
nania izotermy będzie celem następnej części pracy. 
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WYNIKI NUMERYCZNE I DYSKUSJA 


Równanie (17) jest funkcją 0s— 0s (go) dana w postaci uwikłanej. Z tego 
też powodu analiza przebiegu tej funkcji zwykłymi metodami badania 
funkcji byłaby trudna i nieefektywna. Najrozsądniejszą rzeczą jest roz- 
wiązanie tego równania metodą regula falsi przy użyciu maszyn matema- 
tycznych. Uzyskać tak można w efekcie postać tabelaryczną funkcji Qs, 
dla różnych wartości parametrów charakteryzujących oddziaływania mię- 
dzymolekularne i z powierzchnią adsorbentu. Przy obecnej powszechnoś- 


lg 9o 


Ryc. 3 


ci stosowania maszyn matematycznych procedura tego rodzaju nie powin- 
na być uważana za niepraktyczną. Warto wreszcie dodać, ze różne urza- 
dzenia pomocnicze używane w charakterze przystawek komputerowych 
pozwalają uzyskać gotową graficzną postać funkcji o, lub jej obraz na ek- 
ranie oscyloskopu. 
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Ryc. 4 


Gdyby mówić o jakichś poważniejszych zagadnieniach związanych z 
wykorzystaniem równania (17), to należałoby raczej wspomnieć o nie- 
trywialnym problemie obliczania współczynnika f. Współczynnik f dla 
potencjału wyrażonego wzorem (3) został obliczony dość dawno [17] i wy- 
rażony dla jednego mola adsorptywu posiada następującą postać: 

* vtr) 


: R=ZTON GL 1)- i (18) 


Bardziej skomplikowanie natomiast przedstawia się sprawa obliczania fs, 
przy założonych we wzorach (2) i (3) postaciach potencjału. Najważniejsze 
jest to, że wynik całkowania w równaniu (12) na f, zależy od wzajemne- 
go stosunku parametrów: a, b, o, R. 

Lektura najnowszych prac z ostatnich lat skłania . jednak autorów do 
przyjęcia założenia, że b Ko. Crowell i Krizan [18] badając ad- 
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sorpcję na graficie stwierdzili, ze współczynniki wirialne w rozwinięciu 
faktora ściśliwości dla fazy zaadsorbowanej niewiele polepszają opis przy 
przejściu od modelu dwuwymiarowego adsorpcji do modelu trójwymia- 
rowego. Oni też pierwsi potraktowali odchylenia od planarnego układu 
molekuł zaadsorbowanych jako mały parametr. Ten punkt widzenia kon- 
tynuują najpoważniejsze aktualnie prace w tym zakresie, publikowane 
przez Johnsona i Kleina [19] oraz Wolfe i Samsa [20]. 
Jeśli chodzi o wartości liczbowe parametru o, to zostały one już daw- 
no stabelaryzowane [17] dla wielu substancji, a dla najprostszych substan- 
cji (gazów) mieszczą się one w przedziale od 1078 do 1077 cm. Dalsze ba- 
dania wygodnie będzie prowadzić przy użyciu zredukowanych parame- 
trow: U, e oraz y, zdefiniowanych następująco: 
u=% , ene, y= (19) 
Przy użyciu tych zredukowanych parametrów, przyjęciu założenia, że 


Ryc. 5 
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Ryc. 6 
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b Ko oraz zaniedbaniu wyrazów zawierających wyższe niż jeden potęgi 
y. otrzymuje się następujący wynik na f;: 


Bo =$ 61 n[eY (62-1) - e] (20) 


Natomiast fy wyrazi sie następująco: 


B-$6yN|e'(R-1)-R] ŻW 
Dalej równanie (17) przyjmie postać: 
Pr disi iy (22) 


Rachunki numeryczne prowadzone przez nas odbywaly sie wedlug na- 
stepujacego schematu. Badali$my zmiany funkcji o, wraz ze zmianami jed- 
nego tylko parametru, pozostalym parametrom nadajac najbardziej typo- 
we warto$ci. Takie obliczenia robiono dla kazdego parametru z osobna. 
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Wyniki tych obliczeń przedstawiono na ryc. 1—6. Przy odpowiednich 
krzywych pokazano, dla jakich wartości zmienianego parametru się one 
odnoszą, a prócz tego podawano komplet pozostałych ustalonych parame- 
trów. 

Obecnie dokonamy przeglądu najważniejszych własności izoterm 
otrzymanych z równania (17). Dyskusję rozpoczniemy od wpływu średni- 
cy o na postać izotermy adsorpcji. Najbardziej interesujące jest to, że 
wpływ ten zaznacza się wyraźnie dopiero dla małych molekuł, gdzie o < 
< 1077 cm. Widać to wyraźnie na ryc. 1, gdzie postać funkcji o, nie zmie- 
nia się praktycznie dla o < 1077 cm. Z kolei dla o=1078 cm pewne partie 
izotermy, gdzie o, — 0, mają charakter niefizyczny. Pojawienie się takiego 
wyniku nie powinno jednak budzić zdziwienia, jako że równanie (17) po- 
wstało z częściowych sumowań do nieskończoności i ma charakter przy- 
bliżony. W obszarze średnich stężeń oy izotermy zachowują się bardzo ty- 
powo (dla stężeń powyżej 1072 mol/l). 

Układ wzajemny izoterm na ryc. 2 można wyjaśnić dość prosto. Wzrost 
parametru R powoduje wzrost międzycząsteczkowych oddziaływań przy- 
ciągających w warstwie zaadsorbowanej, co ułatwia adsorpcję. Podobny 
efekt powinno dać zwiększenie U, co potwierdza ryc. 3. 

Natomiast trudno wyjaśnić w ten sposób ryc. 4. Warto tylko podkreś- 
lié wyższość izotermy przedstawionej jako funkcja ọs= Qs (oo) nad izoter- 
mą przedstawianą w postaci Na=Na (ge). Funkcje Na=Na (Qo) miałyby w 
tym przypadku podobny przebieg i wzajemne ułożenie. Trudno byłoby 
jednak stwierdzić, czy polepszenie adsorpcji ze wzrostem y wiąże się z pro- 
stym efektem wzrostu objętości, gdzie zachodzi adsorpcja, czy też ma bar- 
dziej skomplikowany mechanizm. Izoterma w postaci os=es (go) daje jed- 
noznaczną odpowiedź na korzyść drugiej z wymienionych możliwości. 

Porównajmy wreszcie ryc. 5 i 6. Wynika z nich, że dla bardzo dużych 
molekuł adsorpcja nie zależy praktycznie od reszty parametrów charak- 
teryzujacych oddziaływania i dla wszystkich tego typu substancji powin- 
na przebiegać w sposób prawie jednakowy. Pozostaje to w zgodzie z wnio- 
skami, jakie zostały wyciągnięte przy okazji studiowania ryc. 1. Są to 
wnioski bardzo intrygujące. Molekułami, dla których o > 1077, powinny 
być na przykład parafiny z liczbą atomów węgla powyżej dziesięciu. Nie- 
stety autorzy nie znaleźli w literaturze odpowiednich danych doświad- 
czalnych dla celów ewentualnej konfrontacji. 
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PE3IOME 


TIpoBoĄuTca uacr4uHoe CYMMHPOBAHME CTaHĄapTHOK pa3BepTKM no cTe- 
NEHN AKTHBHOCTH M30TepMb! AĄCOPÓLMH n no GECKOHEUHOCTH. C 3104 uenbio 
npumeHaeTca rpacbuueckoe uaoGpaxeHue Ypcenna-Mayepa Ana kosdodoa- 
LIMEHTOB pa3aBepTKM, A TAKXKE HEKOTOPbie TeopeTHUECKME pe3ynbTarbl, Nony- 
ueHHbie QynuHbckuM. [lony4eHHbie Takum cnoco60M ypaBHeHMs MW30TepMbl 
aącopóunu oóna4atoT MHOTHMM MHTepecHbIMAM CBOŃCTBAMK, OCOÓEHHO B cny- 
qae Hu3kux HamieHuń aącopGaTa u 6onbiunx Monekyn. 


SUMMARY 


A partial summation of the standard expansion in the activity powers 
of the adsorption isotherm is performed through all terms up to imfinity. 
For this purpose the Mayer-Ursell graphical representation is used, and 
some theoretical results obtained earlier by Fuliński. The isotherm equa- 
tion obtained in this way exhibitis very interesting properties, in parti- 
cular at low densities of adsorbate, and in the case of big molecules. 
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